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1 Introduction 

In recent years, there has been an increasing interest in studying constrained variational problems 
with a fractional diffusion. One of the motivations comes from mathematical finance: jump- 
diffusion processes where incorporated by Merton |14) into the theory of option evaluation to 
introduce discontinuous paths in the dynamics of the stock's prices, in contrast with the classical 
lognormal diffusion model of Black and Scholes [2]. These models allow to take into account large 
price changes, and they have become increasingly popular for modeling market fluctuations, both 
for risk management and option pricing purposes. 

Let us recall that an American option gives its holder the right to buy a stock at a given 
price prior (but not later) than a given time T > 0. If v{t,x) represents the rational price of an 
American option with a payoff at time T > 0, then v will solve (in the viscosity sense) the 
following obstacle problem: 

mm{Cv, V — ^p} = 0, 
v{T) = i;. 

Here Cv is a (backward) parabolic integro-differential operator of the form 



n 1 

Cv = -Vr -rv + ^(r - di)xiVxi - X] ^i^j^v'^^i^j 



2 

i=l itj=^ 



^;(r,a;le^^ . . . ,x„e^") - v{t,x) - ^(e^^ - l)xiVxi{T,x) ^i{dy) 

i=l 



where r > 0, G M, cr = (ujj) is a non- negative definite matrix, and is a jump measure. (We 
refer to the book [S] for an explanation of these models and more references.) When the matrix 
a is uniformly elliptic, after the change of variable Xj i— ?> log(xi) the equation becomes uniformly 
parabolic (backward in time) and the diffusion part dominates. In particular, if no jump part is 
present (i.e., = 0), then the regularity theory is pretty well-understood (see, for instance, |12j). 
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Here we assume that there is no diffusion (i.e., cr = 0), so ah the regularity should come from 
the jump part. We also assume that the jump part behaves, at least at the leading order, as a 
fractional power of the Laplacian, so that the equation takes the form 

Cv = -v^-rv-b-Vu+{-Ayv + K:v, s e (0,1), (1-1) 

where b = {di — r, . . . ,dn — r), and ICv is a non-local operator of lower order with respect to 
{-AYv. 

We now observe that the choice of s S (0, 1) plays a key role: 

- s > 1/2: In this case {—Ayv is the leading term, so the regularity theory for solutions to 
(jl.ip is expected to be the same one as that for the equation 

r min{-i;^ + (-A)^i;, v - ip} = on [0, T] x M", 

\ v{T) =^ onW. ^ ' 

- s < 1/2: If s < 1/2 then the leading term becomes b ■ Vv, and we do not expect to have 
a regularity theory for (jl.ip . On the other hand, in the borderline case s = 1/2 one may 
expect some regularity due to the interplay between b-Wv and — (— A)*-!; (but this becomes 
a very delicate issue). However, when 6 = 0, even if the diffusion term is of lower order 
with respect to the time derivative, the equation is still parabolic and one may hope to 
prove some regularity for all values of s. 

The goal of this paper is to investigate the regularity theory for the model equation ()1.2p . The 
reason for this is three-fold: first of all, considering this model case allows to avoid technicalities 
which may obscure the main ideas behind the regularity theory that we will develop. Moreover, 
since there is no transport term inside the equation, we are able to prove that solutions are as 
smooth as in the elliptic case [U] for all values of s G (0,1). Hence, although when s < 1/2 
the time derivative is of higher order with respect to the elliptic part {—AYv, the regularity of 
solutions is as good as in the stationary case. Finally, as described in Section [5l once the general 
regularity theory for solutions of (|1.2p is established, the adaptation of these proofs to the more 
general case (jl.ip when s > 1/2 should not present any major difficulty. 

Let us remark that the fact that the smoothness of solutions of p.2p is the same as in the 
elliptic case may look surprising. Indeed, the optimal regularity for the stationary problem 
min{(— A)*u, v — ip} = is C^+*(M") [H [151 E|- On the other hand, as we will show in Remark 
13.71 for any /3 G (0, 1) one can find a traveling wave solution to the parabolic obstacle problem 
mm{—VT- + (— A)^/^f,z; — ^} = which is C^~^^ both in space and time, but not C^^"' for any 
7 > /3. Hence, in order to prove that solutions to (|1.2p are in space, one has to exploit the 

crucial fact that v coincides with the obstacle at time T. 

2 Description of the results and structure of the paper 

In this section we introduce more in detail the problem, and describe our main result. 
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Let us observe that, by performing the change of variable t = T — t, all equations introduced 
in the previous section become forward in time. From now on, we will always work with t in 
place of T, so the payoff ■0 becomes the initial condition at time 0. 



2.1 Preliminary definitions 

The fractional Laplacian can be defined as 

i-Arf:=\^f V/gC,°°(]R"), 

so that J f{—A)^g = {f,g)ija- There are also two other different ways to define (—A)*. The 
first one is through an integral kernel: there exists a positive constant Cn s such that 



where the integral has to be intended in the principal value sense. (This can be proved, for 
instance, by computing the Fourier transform of The second one is through a Dirichlet- 

to-Neumann operator, as shown in |7|: given a G (—1,1), for any function / € C^(]R") denote 
by F : M" x M+ M the La-harmonic extension of /, i.e., 

LaF{x, y) := div,,j,(y'^V,,j,F(x, y)) = on x M+, 
F(x,0) = /(x) onR". 

Then there exists a positive constant c„,s such that 

lim y'^Fj,(x,y) = -c„,,(-A)7(x, 0), s := ^ G (0,1). 

In the sequel, we will make use of all of the three above characterizations of the fractional 
Laplacian. However, in order to simplify the notation, we will conventionally assume that Cn,s = 
Cn,s = 1, so that 

-i-^rf= [ TT^^4^dx'= ]im y^Fy{x,y). (2.1) 
J \x' — xl"'^^* y-^o+ 

We will also need the notion of semiconvex function: a function w : ^ R is said to be 
C-semiconvex for some constant C € M if w + C|xp/2 is convex. 

Finally, to measure the regularity of the solutions we will use space-time Holder, Lipschitz, 
and logLipschitz spaces: given a, (3,^,6 € (0, 1), and [a,b] C M, we say that: 
weC^f{[a,b]xR^) if 

11 11 \w(t, x) — wit' , x')\ 

= lkllL-(MxR«) + sup I 1^ < +°°; 
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w G UptCSi[a,b] X M") if 

\w{t, x) — w{t' , x')\ 

J 

[a,b\ X 

w G logLiptCf([a,6] X M") if 



|lf (t, x) — w{t' ^ x')\ 

ll^lllogLip,cf (MlxM^^) ^= lkllL^([a,.]xM«) + ^JUP |t _ + |log |t - t'| | ) + |x - x'l/' ^ 

We will also use the notation w € C"^°'^''^([a, 6] x M") if 

w eCl-''\[a,h]xW) Ve>0, 

and w G ((a,6] x R") if 

G ([a + e,6] x M") Ve > 
(analogous definitions hold for the other spaces). 

2.2 The main result 

Let ^ : M" — )■ M"*" be a globally Lipschitz function of class satisfying (]^_|_||!|'-)L+2s < +00 
and (-A)> G L°°(M"). Fix s G (0, 1), and let u : [0, T] x M" ^ M be a (continuous) viscosity 
solution to the obstacle problem 

min{nt + (-A)^'u, u - ^/;} = on [0, T] x M", 

'u(O) = V on M". ■ 

Existence and uniqueness of such a solution follows by standard results on obstacle problem^. 
The main goal of this paper is to investigate the smoothness of solutions to the above equations, 
planning to address in a future work the regularity of the free boundary. 
Our main result is the following: 

Theorem 2.1. Assume that i) (^C'^{W), with 
IIVVII L°°(IR") + II^Vll 

and let u he the unique continuous viscosity solution of (|2.2p . Then u is globally Lipschitz in 
space-time on [0, T] x M", and satisfies 

r ntGlogLip,Ci-''((0,r] xM"), (-A)^n G logLip,Ci-^((0, T] x M«) if s < 1/3; 
\ ut£ C7^y"°^'^"'((0,T] X M"), {-AYu G C7^y'^^'((0,r] x M") if s > 1/3. 



^ Here, existence of solutions is not the main issue; for instance, one can construct solutions by using proba- 
bilistic formulas involving stochastic processes and stopping times [9]. Another possibility is to approximate the 
equation using a penalization method (as done in the proof of Lemma l3.ip and then use the a priori bounds on 
the approximate solutions (see the proofs of Lemmas 13.21 and I3.3|l to show existence by compactness. The fact 
that these two notions of solutions (the probabilistic one and the one constructed by approximation) coincide, 
follows from standard comparisons principle for viscosity solutions. 
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Let us make some comments. First of all we recall that, for the stationary version of the 
obstacle problem, solutions belong to C^'^'^iW^) (or equivalently, (— A)*ti € C^~^(M'")), and such 
a regularity result is optimal |15| [6], Hence, at least concerning the spatial regularity, our result 
is optimal, too. 

Once the C^~*-regularity of (— A)*it is established, the fact that s = 1/3 plays a special 
role is not surprising: indeed, the operator dt + (—A)* is invariant under the scaling (t, x) i— >■ 

1-3 

(A^*t,Ax). Hence, a spatial regularity Cj:~* naturally corresponds to a time regularity C^^" , 
provided ^ < 1, that is, s > 1/3 (see (|A.3P - (|A.4p in the Appendix). 

Finally, concerning the regularity in time, when s = 1/2 one can construct traveling wave 
solutions which are both in space and time, see Remark 13.71 Hence our result is almost 

optimal in time, at least when s = 1/2 (the result would be optimal if we did not have the 0"^ in 
the Holder exponent). Moreover, the regularity in time is almost optimal also in the limit s — )■ 1 
(since, when s = 1, it is well-known that solutions are in time and C^'^ in space [3l HI |5]). 
Hence, it may be expected that our result is almost optimal in time for all s S (0, 1) (or at least 
for s > 1/3). 

2.3 Structure of the paper 

The paper is structured as follows: first, in Section|3]we discuss some basic properties of solutions 
of (|2.2p . like the validity of a comparison principle, the Lipschitz regularity in space-time, the 
semiconvexity in space, and the boundedness of ( — A)*n. Moreover, we will show that solutions 
are for s > 1/2, and, as explained in Remark 13. 7| C^-regularity in space is optimal when 
s = 1/2 unless one exploits the additional information that the solution coincides with the 
obstacle at the initial time. 

In Section [H we first use an iteration method to show that, for any t > 0, {—Ayu{t) is C" 
near any free boundary point (Subsection 14. ip . Then, we prove a monotonicity formula which 
allows to show that {—A)'^u{t) is C^"'* near any free boundary point for all t > (Subsection 
I4.2p . Finally, combining the fact that ( — A)'^n(t) is C^~** on the contact set with equation ()2.2p . 
a bootstrap argument allows to prove Theorem 12.11 (Subsection I4.3p . 

In Section[5]we briefly describe what are the main modifications to perform in order to extend 
the regularity result in Theorem 12.11 to solutions of (jl.ip when s > 1/2, leaving the details to 
some future work. 

Finally, in the appendix we collect some regularity properties of the fractional heat operator 
dt + {-Ay. 



3 Basic properties of solutions 

Here we discuss some elementary properties of solutions of (|2.2p . Actually, since many of them 
do not rely on the fact that u coincides with the obstacle at time 0, we consider solutions to 

r min{nt + (-A)^'u,'u-^} = on[0,r]xR", 

\ u{0) = no on M", ^ 
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where uq > tp is a globally Lipschitz semiconvex function. Most of the properties of u will be a 
consequence of the following general comparison principle: 

Lemma 3.1 (Comparison principle). Let ipjip : M" — )■ R &e two continuous functions, and 
assume that u, u : [0, T] x M" M are viscosity solutions of 



mm{ut + i-AYu, u - ip} = on [0, T] x M", 
n(0) = no on M", 

and 

mm{ut + (- A)*n, u - ijj} = on [0, T] x M", 
n(0) = no on M", 

respectively. Assume that uq < uq and ip < ip. Then u{t) < u{t) for all t € [0, T]. 



(3.2) 



(3.3) 



Proof. We use a penalization method: it is well-known that solutions of (|3.ip can be constructed 
as a limit of n"^ as e — ?> 0, where is smooth solutions of 

ul + (-A)^n^ = /3,(n^ - ^,) on [0, T] x 

7x^(0) = ng > ^/^e on R", , ^ ' ' 

with n5,V'e e C;?°(R"), l3,{s) = e-^/^ ^ V', (-A)'^^^ ^ (-A)"^/;, and ng ^ no locally 
uniformly as e — ?> (see for instance [HI Chapter 3] for a proof in the classical parabolic case). 

Hence, it suffices to prove the comparison principle at the level of the approximate equations, 
assuming u^{0) < n^(0) and ip^ < ip;,. Let us observe that, since ipi, < and /3'^ < 0, we have 

which implies 

nf + (-A)"n" = /3e(n^ - Ve) on [0, T] x R" 

nf + (-A)^n^ = /3,(n^ - 4) > ^,{u' - on [0, T] x R". 

Since n'^(O) < n^(0), by standard comparison principle for parabolic equations (see for instance 
the argument in the proof of Lemma [3.31 below) we get < u^, as desired. □ 

The following important properties are an immediate consequence of the above result: 



Lemma 3.2. Let u be a solution of (j3.ip . and assume that uq and ip are globally Lipschitz and 
Cq- semiconvex. Then: 

(i) u{t) is Lipschitz for all t ^ [0, T], uiii/i ||Vn(t)||/,oo(i^n) < max{||Vno||j;^oo(-i[j,j), |[VV'||i:,oo(Rn)}• 
('^^j u{t) is Co-semiconvex for all t G [0,T]. 
Moreover, if uq = ip then 

(Hi) [0, r] 9 t I— 7> u{t,x) is non- decreasing in time. 
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Proof, (i) Observe that, for every vector v € M" and any constant C € R, u{t,x + v) + C\v\ 
solves (j3.2p starting from U(){x + v) + C\v\ with obstacle ^/^(x + v) + C|u|. Moreover, if C : = 
max{||Vno||Loo(]Kn), || VV'||l°°{R")}) then uo{x + v) + C\v\ > uo{x) and 7p{x + v) + C\v\ > ijj{x). 
Hence, by Lemma l3.ll we obtain 

u{t, x + v) + C\v\ > u{t, x) Vx, t; e R", t > 0. 

The Lipschitz regularity of u{t) follows. 

(ii) As above, we just remark that u(t, x + v) + u{t, x — v) -\- C\v\'^ solves (|3.2p for every C € R. 
Hence, by choosing C := 2Cq we get Uo{x -\- v) + Uo{x — v) + 2Co\v\'^ > 2uo{x) and ipi^x + v) + 
ip{x -v) + 2Co\v\^ > 2'ip{x), and we conclude as above using Lemma 13.11 

(iii) We observe that, for any e > 0, the function u{t-\-e,x) solves ()2.2p starting from u{e,-). 
Hence, since u(e, •) > ip, by the comparison principle we obtain 

u{t + e,x) >u{t,x) Vt,e>0. 

□ 

We now prove the following important bounds: 
Lemma 3.3. Let u be a solution of ()3.ip . Then 

0<ut + {-AYu < ||(-A)>||ioo(Kn), (3.5) 

ll^t||L°°([0,T]xR") < II(-^)'*^^o||l°°(R"), (3.6) 
In particular {—AYu is bounded, with 

II(-^)''^^IIl=°([0,T]xR") < ll(-A)>||Loo(Kn) + ||(-A)^Uo||Loo(]gn). (3.7) 

Proof. As in the proof of Lemma [3. 1| we use a penalization method: we consider solutions u'^ to 
p.4p . and we prove a uniform (with respect to e) L°°-bound on both (3e{u^ — V'e) ^"^^ uf. 

• L°° -bound on f3^{u^ — ipe)- Since /Sg > 0, we only need an upper bound. 

Assume that inf[Q 'p]xR"(^^~^£) < (otherwise the problem is trivial), and let 93 be a smooth 
function which grows like at infinity. Then, since vanishes at infinity (being a solution to 
a smooth parabolic equation starting from a compactly supported initial datum), for any 5 > 
small we can consider (t^,Xg) a minimum point for — V'e + 7^ "I" ^'f [O'^l ^ -l^"- 

course, min^^^JxR" (^^^ — "^e + + < for (5 sufficiently small, which implies that (i^,xf) 
belongs to the interior of (0, T) x R". Hence 

ul{tl4) + jjA^ = 0, i-AYu^tixi) - {-AYAi4) + 5i-AYv^ixi) < 0, 

which combined with (|3.4p gives 

Piein' - i^sYtlxi) < i-AYM4) - Jrh^ ~ < ||(-A)>,||i^(Kn) + 0{5). 
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Since {u^ — 'tp£){t^, x^) inf[o,T]xR"(^'' ~ V'e) as 5 ^> and < we obtain 

sup f3eiu' - iPe) = lim/3,(n" - < ||(-A)>,||ioo(Kn), 

[0,T]xR" 

SO that (|3.5p follows letting e ^ 0. 

• L°° -bound on nf . We use the same argument as in Lemma 2.1]: differentiating ([3? 
with respect to t we obtain that u;*" := uf solves 

wl + {-AYw' = p',{u^ - ipe)w' on [0, T] x W 
w^O) = -{-AYu^iO) on R". 

Since < and ||t(;^(0) ||ioo = ||(— A)*nQ||ioo, using a maximum principle argument (as above), 
we infer that 

Ikt IIl°°([0,T]xR") = II'"^^IIl°°([0,T]xR") < ll'"^^(0)||L°o(]JJn) = ||( — A)*ttQ||2,oo(]Kn). 

Letting e — ?> we get ()3.6p . as desired. □ 

The above result together with Lemma l3.2( i) gives the following: 
Corollary 3.4 (Lipschitz regularity in space-time). Let u be a solution of (|3.ip . Then 
ll^t|lL°°([o,T]xR") + I|V^IIl°°([o,t]xR") < max{||Vno||ioo(Kn), || V'(/'||loo(r")} + || (-A)''no||ioo(iRn). 

In the sequel we will also need the following result: 
Lemma 3.5. Let u be a solution of (|2.2p with 

IIVV'IIl-'CR") + ll(-A)''^/'||Loo(Kn) + IIVuollLoo(Kn) + ||(-A)'*'Uo||L=o(Kn) < +00, 

and fix to > 0. Then 

< {-Ayu{to) < |l(-A)>||ioo(Kn) for a.e. x G {u{to) = V'}, (3.8) 
{-Ayu{to) < on {u{to) > iP}. (3.9) 

Proof. Let us recall that, thanks to Corollary 13. 4| u is Lipschitz in time. So, Lemma 13.2^ 111) 
gives tij > a.e. 

Moreover, since Ut = a.e. on the contact set {u = ip}, u satisfies 

Ut + {—Ayu = in {u > ip}, Ut = a.e. on {u = ip}, 

which gives 

Ut + i-AYu = {{-Ayu)x{u=t^} 

both in the almost everywhere sense and in the sense of distribution. (Observe that the above 
formula makes sense since {—AYu is a bounded function, see Lemma l3.3p . Hence, u solves the 
smooth parabolic equation ut + {—AYu = /, with / globally bounded and vanishing inside the 
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open set {u > So, a simple application of Duliamel formula shows that u is smooth inside 
{u > tp}. In particular, this fact combined with the non-negativity of Ut implies 

(-A)^n(to) = -ut{to) < on {u{to) > ip}, 

that is, ([32]). 

We now prove (|3.8p . Since ut = a.e. on the contact set {u = tp}, (|3.5p gives 

< (-A)^n < ||(-A)>||loo for a.e. (t, x) G = V}. (3.10) 

Now, to show that the bound < (— A)*ti(to) < ||(— A)*'(/'||j;^oo(-iRn) holds a.e. on M" for every 
to S [0,7"], we observe that the map 

t ^ u{t) G LL(M") 

is uniformly continuous (this is a consequence of the Lipschitz continuity in time, see Corollary 
I3.4p . which together with the uniform bound ()3.5p implies that the map 

i^(-A)Mi)eLL(M") 

is weakly continuous. Thanks to this fact, we easily deduce the desired estimate. Indeed, fix 
e > 0, A C {u(to) = ^} a bounded Borel set, and test (|3.10p against the function X[to-£,to]XA- 
Since the sets {u{t) = ip} are decreasing in time (see Lemma l3.2( iii)). we have [to — e,to] x ^4 C 
{u = tp}, which together with (|3.10p gives 

0< r [ {-AYu<\\{-Ayi;\\L^\A\e. 

Dividing by e and letting e ^ 0, by the weak-L^ continuity of t i— ?> (— A)*n(t) we deduce 

< / {-AYu{tQ) < ||(-A)>||loo 1^1 y Ac {u{to) = V'} Borel bounded, 
J A 

SO that the desired bound follows. □ 

We now show that the uniform semiconvexity of u{t), together with the L°°-bound on 
(—A)^u{t), implies that solutions are in space when s > 1/2 (actually, when s > 1/2, 
by elliptic regularity theory the boundedness of (— A)*«(t) implies that u{t) € Cj^*^"'^^ (M")). As 
we will show in Remark 13.71 below, unless the contact set shrinks in time, this regularity result 
is optimal for s = 1/2. 

Proposition 3.6 (C^-spatial regularity). Let u be a solution of (|3.ip with s S [1/2, 1). Assume 
that uq and tp are semiconvex, and that ||(— A)**uo||Loo(-]Rn) + ||(— A)**^||^cx;(]R7i-) < +00. Then 
u{t) G C^(M") for all t G [0,T]. Moreover the modulus of continuity ofVu depends only on s, 
the semiconvexity constant of uq andip, and on ||(— A)*ito||Loo(]Kn) + ||(— A)'''(/'|[j;^oo(ign). 
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Proof. First of all, we claim that, for every fixed t S [0, T], the map x i— ?> —{—Ayu{t,x) is lower 
semicontinuous. Indeed, recall that if Cq denotes a semiconvexity constant for both uq and tp, 
then u{t) is Co-semiconvex for all t G [0,T] (see Lemma [3.2( 11)). Hence (— A)'*n(t,x) is pointwise 
defined at every x G M", and is given by (see (|2.ip ) 



, , , u(t,x + h) + u(t,x — h) — 2u(t,x) f u(t,x + h) — u(t,x) , 

-AYu(t,x)= / ^ 771—7^ — -^—^dh+ / — — — -^dy 

/ \ ^ / I '}\h\n+2s I , \U\n+2s ^ 

n(t, X + /i) + u{t, x-h) - 2u{t, x) + 2Co\h\'^ 

Bi 2|/i|"+2* 

/■ n(t, a; + /i) - K(t, x) 
-CoC{n,s)+ -—2^ dy 

where C{n,s) := |/ip"""2* (i/i. The last integral in the right hand side is continuous as a 
function of x (since u is continuous). Moreover, since the function inside the first integral is 
continuous in x and non- negative (by the Co-semiconvexity) , the first integral is lower semicon- 
tinuous as a function of x by Fatou's lemma. This proves the claim. 

Now, we remark that — (— A)*ti(t, xq) = +oo whenever xq is a point such that the subdiffer- 
ential of u(t) at xq is not single valued. Indeed, suppose that 

r C 

u{t,x) > ^xo,puP2{x) '■= u{t,xo) + max{pi • {x - xo),p2 ■ {x - xq)} - —\x - xof XBi(xo)(^)> 

for some pi 7^ P2- Then, it is easy to check by a simple explicit computation that 

-(-A)Vxo,pi,p2(2;o) = +00 Vs > 1/2. 
Hence, since u{t) > (pxo,pi,p2 with equality at xq, we get 

{-AYu{t,xo) > (-A)Vxo,pi,P2(^o) = +00. 

However, since — (— A)*u(t) is bounded by ||(— A)'^mo||l°°(ir") + ll(~A)^'0||j;^oo(]Kn-) (see ()3.7p ) and 
it is a lower semicontinuous function, the above inequality is impossible. Thus the subdifferential 
of u(t) at a; is a singleton at every point, i.e., u{t) is C^. Finally, the last part of the statement 
follows by a simple compactness argument. □ 

Remark 3.7. The spatial C^-regularity proved in the above proposition is optimal for s = 1/2. 
Indeed, consider the case n = 1 and = 0, and use the interpretation of the (l/2)-fractional 
Laplacian as the Dirichlet-to-Neumann operator for the harmonic extension, as explained in 
Subsection 12.11 (observe that La = Ax^y when s = 1/2). Then, we look for solutions to the 
problem 

niin{nt -Uy,u} = on [0, T] x M, I'S m 

Ax,yu{t) = on [0,T] X M X M+. ^' 

Let us try to find traveling waves solutions to the above equation, i.e., solutions of the form 
u{t, X, y) = w{at + x,y), with a S M. In this case ut = aUx, so t(;(x, y) has to solve 

awxix,0) - Wy{x,0) = when {?i;(x,0) > 0}, Ci ^2) 

Ax,yW = on M X R+. ^ ' ^ 
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By using the complex variable z = x+iy it is easy to construct C solutions to the above equation: 
if we denote p = \z\ = y/x^+y^ and 6 = arg(z), then wp{x,y) := —p^~^^ sin{{l + /3)6) = 
—lm.{z^~^^) is harmonic in the half-space y > and solves 

wp{x,0) =0 on {x > 0}, 

M^^ 0) > 0, = {wp)y on {x X 0}. 

Observe that Wj^ is of class C^~^^ both in space and time (but not more), and solves (|3.12p 
with a = l/tan(/37r). Since f3 € (0,1) is arbitrary, we cannot expect to prove any uniform 
C^^"-regularity for solutions to (|3.1ip . Thus, the C^-regularity proved in Proposition 13.61 is 
optimal. 

On the other hand, we observe that the case > (i.e., the contact set shrinks in time) 
corresponds to a < 0, or equivalently to /3 > 1/2. Hence, in this case the solutions constructed 
above are at least C]^^^"^ , which is the optimal regularity result for the stationary case [Tl[6]. As 
we will show in the next section, solutions to (|3.ip satisfying > are of class C'^'^'^/'^ in space. 
In particular, by Lemma [3.21 this result applies to solutions of ()2.2p . 



4 Proof of Theorem 12.1 

The strategy of the proof is the following: first in Subsection 14.11 we prove a general C"^"'"^*- 
regularity result in space which, roughly speaking, says the following: let v : — ?> M be a 
semiconvex function which touches from above an obstacle il^ : R'^ — )> M of class C^. Assume 
that (— A)*u is non-positive outside the contact set and non-negative on the contact set. Then 
V detaches from in a C"^"*"^* fashion, for some a = a{s) > universal. In particular, as shown 
in Corollary this implies that {- l^f vx{v=i^} G C^(M"). 

Then, in Subsection 14.21 we use a monotonicity formula to prove the optimal regularity in 
space 

(-A)^t;X{„=^}GCi-^(M"). 

Finally, in Subsection 14.31 we apply the above estimate to any time slice u{t) to prove that 
{— /S)^ u{t)Xu(t)=ii) £ uniformly in time. Then, exploiting (|2.2p and a bootstrap argu- 

ment, we get ()2.3p . 



4.1 A general C"+^''-regularity result. 

In order to underline what are the key elements in the proof, in this and in the next subsection 
we forget about equation ()2.2p . and we work in the foUowinggeneral setting: let v^if) : R" — >• M 
be two globally Lipschitz functions with v > t/j. Assume thao 

(Al) ||L>2V;||i^(Kn) =: Co < +oo; 

^ The smoothness assumption on v inside the open set {v > ip} (see (A5)) is not essential for the proof of the 
regularity of v at free boundary points, but it is only used to avoid some minor technical issues. Anyhow this 
makes no differences for our purposes, since all the following results will be applied to v — u{t) with t > Q, and u 
is smooth inside the open set {u > ip} (see the proof of Lemma |3.5[) . 
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(A2) ||(-A)^Vllc;i-(R") <+oo; 
(A3) v-ip, {-AYv G L°°(R"); 
(A4) V is Co-semiconvex; 

(A5) V is smooth and {—AYv < inside the open set {v > ip}; 
(A6) ||(-A)^V||l-(R") > (-A)^^; > a.e. on {v = t/;}. 

Under these assumptions, we want to show that v is (7"+^* at every free boundary point, with a 
uniform bound. More precisely, we want to prove: 

Theorem 4.1. Let v be as above. Then there exist C > and a G (0, 1), depending on Cq, 
|l(-A)''Vllci-={Rn), llv - V'||l<^{R"), and ||(-A)''i;||ioo(Rn) only, such that 

sup It; - VI < Cr"+2'*, sup |(-A)"i;x{„=^} | < C Vr < 1 (4.1) 

Br{x) Br{x) 

for every x G d{u = ■(/'}. 

Before proving the above result, let us show how it implies the following: 

Corollary 4.2. Let v be as above. Then there exist C' > and q G (0, 1 — s], depending on Cq, 
11^ ~ V'IIl°°(R")> and ||(— A)''u||^oo(iRn) only, such that 

||(-A)''i;x{„=^}||c^(Rn) < C". 

Proof. Without loss of generality, we can assume that the exponent a provided by Theorem 14.11 
is not greater than 1 — s. Moreover, since {—A)^v is bounded on {v = ip} (see (A6)), it suffices 
to control |(— A)**f(xi) — {—Ayv{x2)\ when xi,X2 (z {v = ip} and \xi — X2\ < 1/4. 

Let M := \\v — V'||L°o(Rn)- Moreover, given x £ {v = ijj}, let dpix) denote its distance from 
the free boundary d{v = ip}. 

Fix xi,X2 G = ip}, with \xi — X2I < 1/4. Two cases arise. 
• Case 1: maxj=i^2 dF{xi) > 4|xi — X2I. Set v := v — ip. Since a < 1 — s by assumption, thanks 
to (A2) it suffices to estimate {—AYv inside {v = 0} = {v = ip}. Now, by Theorem 14. II we have 

sup \v\ < Cr''+^', Vr < 1, i = 1,2. 

Br{Xi) 

Hence, since v = inside i?4|^^_^2|(^i) ^ -^4|xi-a'2| (^2) and l^l < M outside Bi{xi) D Bi^2{^2), 
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we get 

\{-Arvix,)-i-Ayv{x2)\ 

V{x') - V{xi) , f V{x') - V{X2) , , 



\x' — Xi 



\n+2s 



dx 



\x' — X2|"~''^'^ 



< 



'R"\[B4|^^_^2|(^l)"^4|a;i-a:2|(a^2)] 

<C\xi- X2I 



\v{x')\ 



dx' 



1 



1 



J\a+2s 



Bi (xi)\[B4|^^_^2 I (^1)^-8412:1 -2:2 I (^2)] 

1 



1 



dx' 



+ 



dx' 



+ 2M\xi-X2\ I 



< C 



c 



\Bi(xi) 

ds + M 



+ 



1 



Xl-X2\ 



Xl — X2\ < C Ixi — X2I", 



dx' 



as desired. 

• Case 2: maxj=i^2 c^Fl^^i) < — X2|- For every i = 1,2, let € = ^} denote a point 
such that \xi — Xi\ = dpixi). Then, by Theorem 14. II we get 

\{-Ayv{x,) - {-Ayv{x2)\ < sup \{-Ayv\+ sup \{-Ayv\ 



B. 



< 8C\xi -X2\ 



□ 



4.1.1 Proof of Theorem 14.11 

The strategy of the proof is analogous to the one used in [1] to study the stationary fractional 
obstacle problem with s = 1/2 (also called "Signorini problem"). 

With no loss of generality, we can assume that is a free boundary point, and we prove (j4.ip 
at a; = 0. Moreover, by a slight abuse of notation, let us still denote by v : M" x — > M the 
La-harmonic extension of v, i.e., 

Lav{x,y) = divx,y{y°'^x,yvix,y)) = for y > 0, 

and v{x,0) = v{x), with v{x) as above. Then 

lim y^Vy^x.y) = -{-Ayv{x,Q), a = 1 - 2s 

(see (|2.ip ). Let us observe that the Co-semiconvexity of u(x,0) (see (A4)) propagates in y: since 

v(x + /i,0) + u(x-/i,0) -2?;(x,0) > -2Co|/iP V/iSM", 
the maximum principle implies 

v{x + h,y) + v{x -h,y)- 2v{x,y) > -2Co|/ip V/i e M", y > 0, 

that is, 
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(A7) v{-,y) is Co-semiconvex for all y > 0. 
In particular, since LaV = we get 
(A8) dyiy^vy) < nCoy". 

In the sequel, we will informally call the above property "a-semiconcavity" in ?J1. Set now 

v{x,y) := v{x,y) - tpix), 

and denote by A := {v{x,0) = 0} = {v{x,0) = ip{x)} the contact set. Observe that Vy = Vy, 
which together with thanks to (A1)-(A6) gives that the function v enjoys the following properties: 

(Bl) v{x, 0) > for aU {x, y) G M" x M+ \ A x {0}. 

(B2) dy{y"'Vy{x,y)) < 2nCo?/", v{-,y) is (2Co)-semiconvex for all y > 0. 

(B3) \imy_^Q+ y°-Vy{x, y) < for a.e. x € A, limj^_j.o+ y'^Vy{x^ y) > for a.e. x S M" \ A. 

(B4) i{x, y) - v{x, 0) < y2 foj. x e A. 

(B5) if v{x,y) > h, then v{x,y) > h — C^p^ in the half-ball 

HBp{x) := {z G Bp{x) C M" : {V^v{x,y),z- x) > 0}. 

Observe that the proof of (B1)-(B5) is almost immediate, except for (B4) for which a (simple) 
computation is needed: using (B2) and (B3), for a.e. x G A we have 

JO Jo ^ Jo s i + a 

and by continuity the above inequality holds for all x G A. 

We use the notation : — x [0, Tj^i a r], where ry^.a = y'^^^' ^^^^ show a decay result 
for y'^Vyi 

Proposition 4.3. There exist two constants Ki > 0, /z G (0, 1), depending on Cq, Iju — V'||l°o(k"), 
and ||(— A)'*'i;||j;^oo(Kn) only, such that 

mfy%y>-Kiti\ (4.2) 
Proof. We prove the result by induction. 

Case k = 1: since v is La-harmonic, y°'Vy = y"'Vy solves the "conjugate" equation L^a{y"'^y) = 
inside R'^xM"*" (see for instance [3 Subsection 2.3]). Hence, the boundedness of limj^_5.o+ y"'Vy{x, y) = 

^ Even if we use the names "a-semiconcavity" and "Co-semiconvexity" with different meanings, this should 
create no confusion. Observe also that, when a = 0, (A8) reduces to the classical notion of semiconcavity. 
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(— A)*f(x,0) (see (A6) and (|2.ip ) combined with the maximum principle impUes the result. 
Induction step: Assume the result is true for k = ko, i.e., 

inf y^vy > -Ki 

for some constants Ki > and € (0, 1) which will be chosen later, and renormalize the solution 
inside Fi by setting 




It will also useful to consider the -harmonic function 




where v is the L^-harmonic function given by 



v{x, y) := v{x, y) — ip{0) — VTp{0) ■ x. 



(4.3) 



Then, thanks to (Al) and (B2): 



(i) \V{x,y)-V{x,y)\ < 



Xl(42(l-»)^)'=0 



X 



^ and Vy = Vyi 



(ii) mir,y''Vy = mir,y%>-l; 



(iii) dyiy'^Vy) = dy{y%) < 




V and V are 



i^i(42( 



2Co 

!{l-s)^)feO 



-semiconvex inside 



Fix L := Cn,aCo, where Cn,a 3> 1 is a large constant depending on n and a only (to be fixed 
later), and define 



W(x,y) := V(x,y) —7- — , , , bp . 



Thanks to (B1)-(B3), the function W satisfies the following properties: 



1. it is L„-harmonic in the interior of F 



i/s; 



2. W{x, 0) < for X € (A \ {0}) x {0}; 

3. lim(^x,y)-^{o,o) W{x,y) = 0; 

4. limy^o+ y^Wyix, y)>0 for x ^ Bi/s \ A. 
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Hence, up to replacing by + ey^ " with e > (so that the inequahty in 4. becomes 
strict) and then letting e — )• 0, by Hopf's Lemma W attains its non-negative maximum on 

dVys \ {y = 0}. 

Two cases arise: 

Case 1: The maximum is attained on dTi/g D {y = T/n,a/8}. 
In this case, there exists xq G -Bi/s such that 

V Vn,a \ ^ _ / ^ 

{^0, 8 ; - '^"'«K^(42(l-s)^)fco' 

for some constant ^ > depending on n, a only. Thanks to the semiconvexity in x (see 
property (iii) above) and recalling that L ^ Cq by assumption, there exists an ra-dimensional 
half-ball HB1/2 {xq, ^) such that 

(see property (B5)). Recall now that limj,_j.o+ y°'^y{Xj v) = li™?/-^-o+ y"'^y{^^ y) > when F(x, 0) > 
0, while \\m.y^Q+ y°'Vy{x,y) = \\m.y^Q+ y'^Vy{x,y) < when V{x,0) = 0. Hence, by the "a- 
semiconcavity" of y in y (property (iii) above) and by (i), it is easy to see that 

^hm y% ix,y) > -C'Lj^^^^^^^yj^ Vx € HB,/, (xo,0) , 

for some universal constant C^'^ > 0. 

Case 2: The maximum is attained on dTi/g \ {y = ??n,a/8}. 

Let (xq,?/q) be a maximum point. Since such a point belongs to the lateral side of the cylinder, 
recalling the definition of r]n,a we have Ixgl^ > i^S^lyol'^i which implies 



Vix'o,y'o)> 



l+a I 

L 



^^(42(i-s)^)fco- 



Again, we recall that limy^Q+ y"'Vy{x,y) = l[my_^Q+ y"'Vy{x,y) > when V{x,0) > 0, while 
limj^_^Q+ y"'Vy{x,y) = limj^__^Q+ y'^Vy{x,y) < when V{x,0) = 0. Thus, by the half-ball estimate 
(B5) applied to V, by the "a-semiconcavity" of y in y (property (iii)) and by (i), we obtain 

lim y% {x,y)>0 Vx € HBy2 {xq, O) . 

Hence, in both case we have reached the following conclusion: 
There exist a constant Ci > 0, depending on n, a, and Cq only, and a point x S -B^/g C M", such 
that 

lun^ y%{x,0) > -^^^^^ Vx E HBy, (x,0) . 
Thus, if we choose Ki and fi satisfying Ki > 2Ci and fi > l/4?^^~''^\ then we obtain 

limy%{x,y)>-l. (4.4) 
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Moreover, thanks to (ii), 

y% > -1 in Fi. (4.5) 

As we already observed before, the fact that V is La-harmonic imphes that y"'Vy solves the 
conjugate equation L_a(y'^V^) = inside x M+. Hence, thanks to (|4.4p and ()4.5p . the Poisson 
representation formula (see [3 Subsection 2.4]) implies the existence of a constant 9 < 1 such 
that 

Therefore, by (i) and the "a-semiconcavity" of u in y (property (iii)), we obtain 
y%{x, y) = ,%(x, y)>-e- j^J^^Z)^^., =■ > "1' 

provided Ki is sufficiently large. Rescaling back, this proves (|4.2p with k = /cq + 1, which 
concludes the proof. □ 

Recalling that y"'Vy = y"'Vy, thanks to ()2.ip and (A6) the above proposition implies 

sup |(-A)^i;X{.=v,}| <Cr° < 1- 

Br{x) 

We now show that a control from below on y"'Vy inside gives a control from both sides on v 
inside T^.i^ . This will conclude the proof of Theorem 14.11 

Lemma 4.4. Fix K > 0, a £ (0,1), and assume that 

infy% >-ii:r" (4.6) 

for some r € (0,1]. Then there exists a constant M = M{K,a,CQ), independent ofr, such that 

sup \v\ < Mr^+"-" = Mr"+2^ 

Proof. Since v is globally bounded, it suffices to prove the result for r small. First of all, let us 
observe that, thanks to (Bl) and (|4.6p . 

v{x,y)>v{x,0)-K r ^du>-K^ V(x,y)Gr„ (4.7) 

Jo u"- 1 - a 

which proves the lower bound on v. 

To prove the upper-bound, assume that there exists a point {x,y) € F^/g such that v{x,y) > 
j^^i+a-a £qj, gQj^g large constant M. Arguing as above, this implies 



provided M is sufficiently large (depending only on K). Now, let B' := B',^„^ar {x, ^"2^ ) denote 

2 

the (n + l)-dimensional ball of radius centered at (x, G M" x M+, and set B'/2 := 

^;„,.r/4 Then B' C F, and (O, ^) € B'/2. 
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Let V be as in ()4.3p . Thanks to (Al), |u — {;| < Cq inside (observe that r]n,a < 1), which 

together with (|4.7p imphes that w + h Cqt^ is non-negative inside Tr- Hence we can 

apply Harnack inequahty inside B' (see [6l Proposition 2.2] and [10]) to obtain 



M 



,1+0— a 



< sup 

B'/2 



A+a—a 



v + K- 



l-a 



< C 



1 - a 



that is t;(0, 7?„^a'"/2) > coAfr^"'"'^ ° — " — Cor^ for some universal constant cq > 0, which 

gives 



v{0,!hh^] >coMr 



l+a—a 



A+a—a 



K- 



2 J - " I- a 

Since € A, combining the above estimate with property (B4) we get 



= i)(0, 0) > V (0, ^) - > coMr 

\ 2 / 1 + a 



which shows that M is universally bounded, as desired. 



„l+a—a 



K- 



1 - a 



1 + a 



□ 



4.2 Towards optimal regularity: a monotonicity formula 

We use the same notation as in the previous subsection. 

We have proved that (—Ayvx{v=4!} grows at most as r" near any free boundary point, 
which implies that {—AyvX{v=tp} ^ C'"(M") (see Corollary 14. 2p . Consider now the function 
w : X M"*" — >■ M obtained by solving the Dirichlet problem 



L-aW = on 
w{x,0) = {-Ayv{x)x{v=i,}{x) on 



(4.9) 



Since w{x,0) > 0, the maximum principle implies w >0 everywhere. 

Assume that G M" is a free boundary point. Since {—AYv{x) is globally bounded (see 
(A3)), using the Poisson representation formula for w [TJ Subsection 2.4] together with the 
uniform C"-regularity of w{x,0) (Corollarv 14. 2p we get 

sup w^x, y) < C r", 

\x\^ -\-y^ <r^ 

for some uniform constant C. The goal of this subsection is to show that 



sup w{x, y) <C ^ 

|xP+y2<f.2 



and 



(4.10) 

-^Yv\\l 



for some constant C > 0, depending on Cq, |1( — A)'''0||^i 
only. 

This estimate will imply that (— A)*t; grows at most as at every free boundary point, 

so that the same proof as in Corollarv 14.21 will give that {—A.Yvx{v=^i)} ^ Ci~*(I^")i with a 
uniform bound. Then, in the next subsection we will apply this estimate to v = u{t) for every 



18 



t G (0,T], and using (|2.2p we will obtain the desired regularity result for u. 



As in the previous subsection, we consider the function v{x,y) = v{x,y) — tjj{x). Thanks to 
Theorem 14.11 together with the (2Co)-semiconvexity of v (see (B2) in the previous subsection), 
we can mimic the proof of [H Lemma 5]: 

Lemma 4.5. Let (7 > and a G (0, 1 — s] be as in Theorem \4-l\ and set 5a = Sa{s) : = 
i (^ a+2s ~ f ) • ^^^^ there exists tq = ro(a, s, C, Cq) > such that the convex hull of the set 
{x G M" : w{x, 0) > r"^"^"} in Br C M" does not contain the origin for r < r^. 

Proof. Thanks to (B3), 

v{x,0)=0 and lim y%(t, x, y) = -'u;(x, 0) < -r"+'^" Vx G {^(x, 0) > r"+'^"}. 

Hence, by the "a-semiconcavity" (B2) of v in y, for any x G {w{x, 0) > r""*""^"} we have 

, X , s'^Wj;(t,x,s) 
v{x,h) < / ^ds 

JO -s" 

1 — a 1 + a 2s 1 + a 

On the other hand. Theorem 14.11 gives 

v{0, h) = v{0, h) - v{0, 0) > -C'/i"+2^ (4.12) 

Assume now by contradiction that the convex hull of the set {(x,0) : w{x, 0) > r"+'^«} n Br 
contains (0,0). Then, by the (2Co)-semiconvexity of v{-,h) (see (B2)) we get 

v{0,h) < sup v{x,h) + CQr^, 

which together with (|4.12p and ()4.1ip gives 

(7/,-+2s > l^a+<5./,2. _ ^^2 _ ^ 2 

- 2s 1 + a ' 

for all r,h £ (0, 1). To get a contradiction from the above inequality, we want to choose h = h[r) 
in such a way that 

/i^ < < < r"+^" /i^'* for r sufficiently small. 

To this aim, set h = r^+^^a/". Then /i" = r"+^''" = o(r"+^"), and both the first and the third 
condition above hold. To ensure that also the second one is satisfied, it suffices to have 

(a + 2s)n + 2^ j < 2, 
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that is 

\ ( a a 



2\a + 2s 2 

Recalling that q + 2s < 2 (so, the right hand side is positive) and the definition of 6a, we get 
the desired contradiction, which concludes the proof. □ 

We now want to use a monotonicity formula to improve the decay of w{x,y) at the origin. 
We first need some preliminary results: 

Lemma 4.6. (i) There exists a constant C' , depending on Cq , |[(— A)*^||^i-s^jg„^, ||^^~V'IIl°°(R"); 
and IK—AyvWi^oof^^n-j only, such that 

(ii) It holds 

Proof, (i) To show the estimate, let us observe that: 

(1) Since 'w{-,0) = on R" \ A while w{-,0) > on A (see (A6)), by the maximum principle 
we get w{x, y) > 0. Hence 

w{x,y) >w{x,0) VxGM"\A, y>0. 

(2) By the a-semiconcavity of in y (see (A8)), 

y^Vyix, y) < lim s^Vyix, s) + 

(Observe that the above limit always exists, since {—Ayv{x,0) is Holder continuous on 
the contact set, while v is smooth outside, see (A5).) 

(3) The function y^Vy solves 

L.aiy^vy) = 0, 

limj^^o+ y''vy{x,0) = -{-Ayv{x,0), 

(see [3 Subsection 2.3]). Since w{x,0) > ( — A)**'y(x, 0) by (A5), the maximum principle 
gives w > —y'^Vy on x Hence, since w{x,0) = —limy^Q+ y"'Vy{x,0) in A, by (2) 
above we get 

w(x, y) > w(x, 0) -y^^" Vx G A, ?/ > 0. 

1 + a 
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Combining (1) and (3) we obtain 



w{x,y) > w{x,0) - -^^^y 
1 + a 



l+a 



Vx € M", y > 0. 



(4.13) 



This estimate, together with the regularity of w and the fact that w is non-negative, implies 
that, for all x ^ and y > 0, 

uP'{x,y) -n?{x,{)) = [w{x,y) - w{x,Q)]\w{x,y) +w{x,Q)] > -Ky^^"-{r + y)"" , (4.14) 

for some uniform constant K > 0. 

We now want to estimate from below 



limsup / {w\{x,y) \un~l^a)/2 
y^o+ J Br (|x|^ + 2/^)1." ^ "-ii^ 



dx. 



To this aim, consider the change of variable s = s{y) := ^'^d define w{x^s{y)) 

w{x,y). Then ()4.14p becomes 



l+a 



r,2 



w 



{x, s) - w'^{x, 0) > -K's{r + si/(^+"))° Vx€Br,s>0, 



(4.15) 



for some uniform constant K' > 0. Moreover, since y '^{w'^)y{x,y) = {w'^)s{x, s), we are left 
with estimating 



limsup / {w'^)s{x, s)- 

S^0+ J Br ' 



dx. 



(|x|2 + (1 + a)2s2/{l+a))(n-l-a)/2 

To do this, we average the above expression with respect to s G [0, e] and we use Fubini Theorem 
to get 



ds I {w )s{x,s) 

J Br 

1 

Br ^ 



dx 



(|x|2 + (1 + a)2s2/(l+»))(«-l"«)/2 

zZ'2(x,e) u;2(x,0) 

-r\n~l~a 



(|x|2 + (1 + a)2e2/(l+a))(n-l-a)/2 



dx 



1 r J /■ ~2, 

- ds w [x,s) — 
e Jo J Br ds 



1 



|2;|2 + (1 + a)2g2/{l+a)){n-l-a)/2 

Now, thanks to ()4.15p . the C"-regularity of w{x,0) = w{x,0), and the fact that 



dx. 



d 



1 



ds V(|x|2 + (1 + a)2s2/(l+a)){n-l-a)/2 

we obtain that the above expression is bounded from below by 
1 



w2(x,0) -K'e(r + ei/(i+'^))" w'^{x,Q) 



\x\'i + (1 + a)2e2/(l+a)){n-l-a)/2 



<0, 



dx 



+ C ' 



1 



Br F 



n— 1— a 



dx 



Br £ 



|x|2 + (1 + a)2e2/(l+a))(n-l-a)/2 



In— 1— a 



dx 
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Concerning the first term in the right hand side, since a = 1 — 2s < 1 we have 

(r + e^/(^+'^) r [ I dx ^ C„,„r"+^+'^ = C„.,r"+i+'^ as e ^ 0. 

For the second term, we want to prove that it converges to as e —t- 0. To this aim, we spht the 
integral into two terms: the integral over B^p, and the one over Br \ B^p, where /? > has to be 
chosen. For the first term, we can bound it from below by 



C 



1 2a 



\n—l—a 



(ix = C7e'^(2a+a+l)~l_ 



Thus, by choosing /3 S y 2a+a+i ' iTaj ensure that the above expression converges to as 
e — 7- 0. Moreover, the fact that /3 < 1/(1 + o) implies that 



Therefore, for estimating the second part we can use polar coordinates and the fact that 

+ (1 + a)2£2/(^+'^))("-l-'^)/2 _ |^|n-l-a^^^2/(l+a)|^|n-3-a y |^| > 

to write 



1 2a 



Br\B p 



e 



1 



(|a;|2 + (1 + a)2e2/(l+a))(n-l-a)/2 



In— 1— a 



dx 



C 

^ — I P 



n-l+2a 



pTi—l—a _j_ (J^2/{l+a) pU—Z—a ^n— 1— a 



dp 



C 

- / P 



2a+a 



p2 + Ce2/(i+«) 



dp = — — 

e 



„2a+a 



e 



2/(l+a) 



p2 + Ce2/(i+«) 



dp 



> 



Q^2/(l+a) rr 



r p2a+a-2 > _^g2/{l+a)-l U ^ ^ 



/3(2a+a-l)l 



Let us remark that 2/(1 + a) > 1, soif2Q; + a — 1 >0 the above expression obviously converges 
to 0. On the other hand, if 2a + a — 1 < 0, since /3 < 1/(1 + a) we get 



^ l + /3(2a + a- 1) > ^ 



2a + a- l 2-l-a + 2a + a- l 
1 H > 



1 + a 1 + a 1 + a 1 + a 

and again the above expression converges to 0. All in all, we conclude that 



2a 
1 + a 



> 0, 



1 2a 



Br ^ 



1 



1 



|2 ^ (1 _)_ Q^2g2/(l+a)^{n-l-a)/2 j^jjn-l-a 



SO that combining all our estimates together we obtain 



1 



£^0 e 



liminf- I ds I {w'^)s{x, s) 



1 



Br 



(|x|2 + (1 + a)2s2/(l+'^))('^-l-a)/2 



dx ^ as e ^ 0, 



dx > -K'C„ar^^^^\ 
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From this fact we easily deduce that 

hmsup / {W%ix,e) ^^^^^ ^ ^ j2/(l+a))(n-l-a)/2 ^ -/^'Cn,ar■"+^+^ 



e->0 J _B. 



which concludes the proof of (i). 

(ii) In this case, we use the C°-regularity of w to control the integral by 

y^~<^ [ I dx 

(|x|2+y2){n+l-a)/2-a " 

Using polar coordinates, the above integral is comparable to 

1 r I f" p"~^ y^~"' 2 

y (p2 _^y2)(n+l-a)/2-Q y + ^^„+i_a-2a ~ yl-a-2a = ' 

and the above expression converges to as — t- 0. □ 

We will also need a result on the first eigenvalue of a weighted Laplacian on the half-sphere. 
We use 'Vqw to denote the derivative of w with respect to the angular variables:. 

Lemma 4.7. Set S^+^ := S'^+i n {x„+i > 0}, := aS!^+^ = S"+i n {x„+i = 0}, := 

S'^+i n {rc„+i = 0} n > 0}. T/ien 

inf = (l-s)(n-l + s). 

ft6/fi/2(§o+^),/i=o on Jgn+i n^y "-da 

Proof. Let denote the restriction to S"^^ of 

H{x,y) := {^/xlTy^ - XnY~% 

that is H = r^~^h{9). As shown in [6l Proposition 5.4], /i is the first eigenfunction corresponding 
to the above minimization problem. Moreover, H solves L-aH = for y > 00 Let Ai denote 
the eigenvalue corresponding to h, so that 

/gn+i \Veh\'^y~"' da 

inf — ^ ro~z — ; — = 

h6Hi/2(§J}+i),h=o on j^n+i h^y "-da 



''a simple way to check this fact is the following: the function G := + — a-„)^''^ is harmonic inside 

y > 0, since it is equal to the imaginary part of the holomorphic function z i— >■ z^^^ , z = Xn + iyn- Moreover, by 
a direct computation it is easily checked that G satisfies 

I V.GI" + {Gyf - ^ =0 y{x,y)£ R" X R+. 

y 



Thanks to this fact, since H = G^'^"") = G^+", we get 

L-aS ^ L.aiG"^^") = (l + a)G''A,,yG + (l + a)aG'' 

as desired. 



|V.Gr + (G,)^-^ 



0, 
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In order to explicitly compute Ai, we observe that divg{y ""X/gh) = Xih. In particular, evaluating 
the above identity at the point (0, 1) G x we obtain 

Ae/i(0,l) = Ai/i(0,l). 

We now write the equation L^aH = in spherical coordinates: 

— Tl — 1 — CL — 

ArH + -Hr + ^AgH - -Hy = 

Evaluating the above expression at (0, 1) € M" X M+ and recahing that H = r^^^h, we get 

= A^i7(0, 1) + nHr{<d, 1) + Ae^(0, 1) - ai7,.(0, 1) 
= -(1 - s)s h{0, 1) + (1 - s){n - a) h{0, 1) + Agh{0, 1). 

Hence 

Agh{0, 1) = -(1 - s)(n - a - s) h{0, 1) = -(1 - s)(n - 1 + s) h{0, 1). 
which gives Ai = — (1 — s)(n — 1 + s) as desired. □ 

To simplify notation, we use the variable z to denote a point {x,y) € M" x IR+. 
Lemma 4.8. Let w and tq be as above, set := {z € M" x R"*" : \z\ < r}, and define 

^(r := — — — / \ y\ dz Vr<l. 

There exists a constant C" , depending onCo, \\{—A)'^7p\\^i~s^^ny 1^— V'IIl°°(R"); '^'^'^ II(~^)'*^IIl°°(R") 
only, such that 



ifir) < C" 



Vr < 1. 



Proof. First of all, we show that <^(1) is universally bounded, so that in particular (^(r) is well- 
defined for all r G (0, 1]. 

Set Lpe{r) := ^317=77/5+ n{y>e} '^2|i-i-£t monotone convergence theorem, it suf- 

fices to estimate liminfe^i^o fei^)- Let x ■ I^" ~^ [0, 1] be a smooth compactly supported function 
such that X = 1 on i?i C M". Then 

Since L-aW = we have L-aiw^) = 2wL^aW + 2\'S/w'^y~'^ = 2\SIw'\^y~°' , which implies that the 
right hand side is equal to 

L liF^^^"^ dxdy = j" j^^ V{w') . V (^^p^) y-x{x) dx dy 

+ [ [ V^(-w^) • V3;x(a;) I ^^_^_^ dxdy 

+ _^ x{x)dx . 



y=e 
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Integrating by parts again the first two terms in the right hand side, and using that L_a jTpr^ 
Cn,aS(ofi), we find that the above expression coincides with 



?/=! r y-a 

+ l_{w\^r—YZ^x{x)dx 



y=e 



y=l 



y=e 



Now, since x = 1 inside the first two terms above are immediately seen to be bounded. 
Concerning the last two terms, the integrals evaluated at y = 1 are clearly finite (and universally 
bounded), since w is smooth for y > 0. Finally, we apply Lemma 14.61 to estimate the integrals 
at y = £, and we obtain 

= liminf < C^, 

£-S>0 

for some constant depending on Cq, ||(— A)'*'0||(7i-s(iRn), lb — V'IIl°°(R")j a-i^d ||(— A)*t)||ioo(]{jn) 
only. Observe that, as a consequence of the fact that (p(l) is finite (i.e., ^"^^^jij^^ — is integrable 
over Bf), we deduce that (/?e(r) — t- ip{r) locally uniformly over (0, 1]. 

Now that we have proved that ip{r) is well-defined, we want to estimates from below its 
derivative. Again, we will do our computations with ip^, and then we let e ^ (jfj. 

Let us assume r > e, and split d {B^ f] {y > e}) as the union of dB^ri{y = e} and dB^r){y > 



^ The proof of the monotonicity formula may look a bit tedious, since we always prove the result at the e level, 
and then we show that one can take the limit as e — > 0. Let us point out that this level of precision is actually 



needed: indeed, assume that we had chosen a different operator Lb (6 G ( — 1, 1)) to define w in (|4.9|l . and we 
defined fir) replacing —a by b (changing, of course, the value of s correspondingly). Then, if one does a "formal" 
proof of the monotonicity formula, one would obtain (at least in the stationary case, so that w{x, 0) = {—A)''v{x)) 
that Lemma [4.8l is true with b in place of —a, and this would imply a false Holder regularity for w (since we know 
that w should be only C^^'^). The fact that we have chosen the "right" operator L-a to define w has played a 
key role in the proof of Lemma 14.61 which is now providing to us some fundamental estimates, which are needed 
to give a rigorous proof of the monotonicity formula. 
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e}. Using again that L^a |^|n-i-a = Cn, 0(^(0,0) ^^'^ recalling that a = 1 — 2s, we easily get 



^l+2(l-s) r" J9B+n{y>e} 

ID Vif • T-; — ^ — da 



1-s r 0. „ / 1 \ , 1 



r 



^ / V(u;^)-vf , , \ ]y-^dz + — [ \Vw{z)\^y-'' da 



2(1-5) /■ ^ ^ 

' —If vw ■ V — — - — da 



l-s 



w'V[j-—^]-i.y--da 



\ [ \Vw{z)\h-''da 
r'' JdB+n{y>e} 

2(1-S) /■ _a , /■ . 2^ 

n — / wWrV da -\ — , , — r / t(; L t-; — ^ — 

(1 - s)(n - 1 - a) /■ 2 -a r 

^ I w^y da 

r"^^^ JdB+n{y>e} 



da 



Thanks to Lemma 14.61 we can estimate from below both the second and the last but one term 
in the last expression. So, letting e — >■ and using that kp^ ^ Lp locally uniformly, we deduce 
that the distributional derivative D^ip of Lp is bounded from below by 



[ wwry-^da -Cr''+^+'' 

[ w^y-da + \f \Vw{z)\''y"'da, 

JdBr.4^ ^ JdBr.^ 



'dBr,+ 

(1 - s)(n - 1 - a) 



for some universal constant C. Now, by Schwartz's inequality the first term in the above expres- 
sion can be estimated from below by 

Hence, recalling that |Vt(;(z)p = {wr)"^ + :pj\V gWgl"^ and observing that n — 1 — a + 1 — s = n — l + s, 
we obtain 

Drip>4T2 [ |V,Mz)|V"da-il^^l%^^i+^ / w'y-^da-Cr-+'+'^. 
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Consider now the function W := {w - r Then \VeW{z)\^ < \Vew{z)\'^. Moreover, by 
Lemma 14.51 ^ is admissible for the eigenvalue problem in Lemma 14.71 Hence 

Dr^>4r2f \^eW{z)\-y~^da-^^^^^=^ j W\-^ da 

r"'^'' JdBr,+ JdBr,+ 



r""^^ JdBr.+ 
> ^^~'^^n+2^^'^ I [(^ - + ^^(^ - ^)] y'^da - C 

^ J dBr + 



Since \W\ < 1^1 < C r"^ and \w — W\ < r""*""^" we obtain 

which integrated over [r, 1] gives 

¥?(r) < +Cr2"+'^"-"-^ +C Vr<L 

(Recall that 1 + a > 0.) Since is universally bounded, this concludes the proof. □ 

We are now ready to prove the optimal decay rate around free boundary points. 

Proposition 4.9. There exists a constant C > 0^ depending on Cq, ||( — ^)^ip\\^i—s ^-^^y \\v — 
'0I!l°°(R"); ^iT-d II ( — ^)'*^^||l°°(M") only, such that ()4.10p holds. 

Proof. Define Ws = w*ps, where ps = Pei^) is a smooth convolution kernel. Since L_a commutes 
with convolution in the x variable, Ws is L_a-harmonic on M" x R'^ . Moreover, by (|4.13p . 
We{x,y) - we{x,0) > -^y^+''. 

Set Ws := (t^e — r""*""^")"*". Then it is easily seen that is L_a-subharmonic for y > 0, and 
Ws{x,y) — We{x,0) > — Y^y^^°. Consider now the function 



Mx, y) ■■= Weix, \y\) + (^1 + ^) \y\^+'^ 



on K X 



We observe that We is L_a-subharmonic outside {y = 0}. Moreover, since Ws{x,y) — Ws{x,0) > 
|y|^+" and Wc is smooth in the x variable, we deduce that it is L„a-subharmonic on the whole 
M" X M. Letting e we obtain that 



w 



(x, y) := {wix, \y\) - r-+'5«)+ + (^1 + ^) |y| 



l+a 



is globally L_a-subharmonic. Thanks to Lemma 14.81 w vanishes on more than half of the n- 
dimensional disc Br x {0}. So we can apply a weighted Poincare inequality (see |10) ) and the 
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definition of ip (see Lemma I4.8p to get 



\Vw\^y'''dz 



J 

< Cr"+2[(^(r) + ri+"] < Cr"+2 [(^(r) + l] Vr < 1, 
which combined with the L_a-subharmonicity of w and Lemma [4.81 gives 

C 



sup(?Z;)^ < 



f.n+1—a 



r/2 



B7 



< Cri+"[(/p(r) + 1] < C 
Hence, since 1 + a = 2(1 — s) we obtain 



j,l+a _|_ ^2a+5a 



sup w < C 

K/2 



sup w + r"+''" + r^+" 



r/2 



< c 



_|_ ^a+5a/2 



Vr < 1. 



Since the above bound holds at every free boundary point, by the very same argument as in the 
proof of Corollary 14.21 we obtain that ||7u||^,8a < C, where /3a = /3a(s) := min{a + Jci/2, (1 — 
s)}. Observe now that, by the formula for 6a provided in Lemma 14.51 given oq > there exists 
6q > such that 5q > (5o > for a S [qq, 1 — s]. Hence, by iterating the above argument k times 
we get 



sup w < C 



B- 



Vr < 1, 



/2 



□ 



and after finitely many iterations we obtain ()4.10p . 

Arguing as in the proof of CoroUarv 14.21 ()4.10p gives: 

Corollary 4.10. There exists a constant C" > 0, depending on Cq, ||(— A)'^'0||^i 
V'I1l°°(ir")> ^iT'd ||(— A)**^^||l°°(R") only, such that 

\\{-AYvx{v=ip}\\ci-''(Rn-^ < C". 

4.3 Almost optimal regularity of solutions to the parabolic fractional obstacle 
problem 

Let ti be a solution of ()2.2p . with ^ G C^(M") satisfying assumptions (A1)-(A2) of Subsection 
14.11 Let us remark that n(0) = ^, while for all t > we can apply the results of the previous 
subsections with v = u{t). Hence Corollary 14.101 gives: 
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Proposition 4.11. Let u^ip he as above. Then there exists a constant Ct > 0, depending on T, 



\D 



2„ 



D(^n), and ||(— A)*V'llci-''(Kn) only, such that 



sup ||(-A)%(t)x{„(t)=^}||ci--(Rn) < Ct. 
fe[o,T] ^ ' 

Proof. As explained at the beginning of Subsection 14.11 u{t) satisfies assumptions (A5)-(A6) 
of Subsection 14.11 for every t > (see Lemma l3.5p . Moreover, (A3)-(A4) follow from Lemma 
13.31 (since \\u{t) — '(/;||^oo(ign) < THuf ||2,oo([o,t]xIR") ^ ^II(~^)''V'IIl°°(R"))- Hence the result is an 
immediate consequence of Corollarv 14 . 1 Ul applied to f = u{t) for any t > 0. (For t = the result 
is trivial since u(0) = □ 

Now, we want to exploit the fact that u solves the parabolic equation 

ut + {-^fu = {{-/\yu)x{u=i,} on (0, T] X M". (4.16) 

Thanks to Proposition HJH the right-hand side of (I44B belongs to L~([0,r]; C^-^(M")), which 
by parabolic regularity implies 

ut,{-^yu e L~((0,T];Ci-^-o^(]R'^)), (4.17) 

see ()A.7p . We now want to use ()4.16p and a bootstrap argument to obtain the desired Holder 
regularity in time. We start with a preliminary result: 

Lemma 4.12. Let that (-A)^^ X{n=v} ^ L°°((0, T]; C^-^(M")). Fix a € [0,min {l, ^}), and 
assume that: 

- Ut e L~((0,r] X M"), {-AYu € L~((0,r];Ci-^-o+(M")) ifa = 0; 

- Ut e C";^""((0,T] X M"), {-AYu € L°°((0,T];C7i-«(M")) ifa>0. 
Then 



(-A)'nx{„=^} G 



C,;:^' ' ((0,r]xR-) ^/a = 0, 

LipiC^-^((0,r] X M*^) «/a>0, s<l/3, 

logLipfCi-"((0,r] X M") ^/a > 0, s = 1/3, 

Cll^"^^'^~'((0,r] X M") tfa>0,s> 1/3, 



a uniform bound. 
Moreover, for s > 1/3 the function 

I — s 

a iH> ^(a) := (1 + a) 

1 + s 

is strictly increasing on [O, , o-nd <I> (^^) = 
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Proof. We need to estimate 

\{-Ayu{t,x)x{u{t)=i,} - {-^yu{s,x)x{uis)=^}\, 0<s<t. (4.18) 

Since {u{t) = ^} C {u{s) = ?/^} for s < t, we can assume that x € {u{s) = "0} (otherwise the 
above expression vanishes and there is nothing to prove). Moreover, since {—AYu vanishes on 
the free boundary, if x S {u{s) = {u{t) = ip} we can alway find a time r € (s, t) such that 

{-Ayu{T,x)x{uiT)=,p} = {-^Yu{t^^)X{u(t)=^p} and x G a{n(r) = ip}. 

Then, if we can estimate (|4.18p with r in place of t, then we will also get the desired bound by 
simply replacing r with t. Hence, we only need to consider the case x G {u{t) = ip}. 

We have to estimate \{—Ayu{t,x) — {—AYu{s,x)\. Let (j) he a, smooth non-negative cut-off 
function supported in Bi such that 0=1, set (/>r(x) := ^0 (^), and compute 



\i-Ayu{t,x) - {-Ayu{s,x)\ < 



+ 



+ 



[{-Ayu{t, x) - {-Ayu{t, z)](l)rix - z) dz 

[{-Ayu{t, z) - {-Ayu{s, z)]<f>r{x - z) dz 
[(-A)*n(s, x) - {-Ayu{s, z)](l)r{x - z) dz 



(4.19) 



We now distinguish between two cases: 

• a = 0. Thanks to the C^~*~'^^-regularity of (— A)*n and the fact that supp (pr C B^, the first 
and the third term in the right hand side of (|4.19p are bounded by C r^~^~^^ . For the second 
term, we integrate (—A)'* by parts, and using that |[( — A)'^0r||Li(R") = II(~^)*0IIli(IR")/'"^'^ ^ 



1/r and that u is Lipschitz in time (Corollary [33]), we get 



\{-Kyu{t,x) - {-Ayu{s,x)\ < c 

Choosing r := \t — we obtain 

\{-Ayu{t,x) - {-Ayu{s,x)\ < C{t - s)TTf-o+, 

as desired. 

• a > 0. Arguing as above, the first and the third term in the right hand side of (|4.19p are 
bounded by C r^~*. For the second one, we integrate again by parts and we estimate 

[{-Ayu{t, z) - {-Ayu{s, z)]4>r{x - z) dz 
[u{t, z) — u{s, z)]{—Ay(j)r{x — z) dz 
\u{t, z) — u{s, z) — ut{s, z)[t — s]\ {—Ay(j)r{x — z) dz 



< 



\ut{s,z)\ j(-A)^0^(x - z)\dz 
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Since ||(— A)'*(/>r||^i(]gn) ~ and ut € C°, the first term in the right hand side is bounded 

by C ^^—^3 — . For the second term, we observe that ut vanishes at {t,x) G {n = -0}, so by the 
C^~'^-regularity of ut we get 



\utis,z)\ li-AYM^- z)\dz 



< C 



mm 



\{-l\ycl)r{x- Z)\dz 



We now remark that, since (j) is compactly supported, \{— l^Y cj){w)\ < ^^yl+2s for \w\ large. So, 
there exists a constant C^, depending on cj) only, such that 

|(-Ar./>H| < V«;GM^ 
1 + 

which by scaling gives 

a 

Hence 



mm 



|(-A)"0^H|du; < C /" ^J^^^, ' ^ dw + C f dw 

I \w\^~'dw + C [ \w\^"''^'-''dw + C, 



j.n+2s 

which implies 



Bl\Br 



mm 



C ifs<l/3; 
{\w\^-%l}\{-AyMw)\dw < { C(l + jlog(r)|) if s = 1/3; 

C{l + r^-^') ifs>l/3. 



\i-Ayu{t,x) - {-Ayu{s,x)\ < c 



All in all, we have obtained 

u _ „^l+a [ C if s < 1/3; 

-'+^- 2^ + {t-s){ C(l + |log(r)|) ifs = l/3; 

^ y C(l + ri"3^) if s > 1/3. 

Choosing r := (t — s)(i+°)/(i+^) , the above estimates give: 
- s < 1/3: \{-AYu{t, x) - (-A)^n(s, x)\ < C{t - s). 

1/3: \{-Ayu{t, x) - i-Ayu{s, x)\ < C{t - s) [l + | log(t - s)\\. 
s > 1/3: Since a < (1 — s)/2s by assumption, we have a < < 1 + ^— — ^q;^-^, so 



- s 



\{-Ayu{t,x) - {-Ayu{s,x)\ < c 



(1 + q)(1-«) , (l-3s)(l + a) 

{t - s) ~' +{t-s) + {t- sY^ ~ 



<C{t-s) — — — , 
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□ 

Thanks to the above lemma, we can use ()4.16p and a bootstrap argument to prove our main 
regularity result. 

Proof of Theorem \2.1\ The global Lipschitz regularity of u in space-time follows from Corollary 

m 

By Proposition 14.111 and (|4.17p . we can apply Lemma [4.121 with a = to deduce that the right 
hand side of (|4.16p belongs to C/^'' ' ((0, T] x M"). Hence, since 2s < l + s, by the parabolic 

regularity theory for dt + (-A)'' (see dTTTll 'l we get ut, {-AYu £ C^~^^'^~' {{0,T] x M"). We 
now apply Lemma |4. 121 with a > 0, and we distinguish between two cases: 

- s < 1/3: In this case we get 

{-Ayux{u=4,} e logLip,Ci-^((0,r] X M"), 

so by (jA.ip and ()A.6p we get {—AYu G logLipjC^~'*((0, T] x M"), and we conclude by using 
Ut = {-Ayux{u=^} - (-A)^^. 

- s > 1/3: Lemma [4.121 gives 

{-Ayux{u=^} G <i^"°^)''"'((o,r] X M"), 

which by ([ATT]) implies ut, {-AYu G c2^~°^'*'^"'((0, T] x M") (recall that $(a) < ^ 
if a < ^^). Hence, we can use iteratively Lemma |4 . 1 2 1 and (jA.ip to get 

ni,(-ArnGC,,,^^+= ^' ((0,r]xM"), 
which together with (jA.SP implies 

(-Art.GQy''-^((0,r]xM"). 



Finally, since ( jj^ — 0+ ) ^^-^ as n ^ oo, we obtain 



i+s ; / 2s 

O+.l-s 



as desired. 



nt^c',^/ " " ^((0,r] xM"), 



□ 
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5 Extension to more general equations 

In this section we give a brief informal description of the main modifications needed to extend 
the regularity result of Theorem 12.11 to solutions of (jl.ip . at least when s > 1/2. Our aim is 
only to point out the major differences with respect to the model case ()2.2p treated above, and 
to explain how to handle them. There will be however to attempt to state a proper theorem, 
as this would need a careful analysis of the assumptions needed on ■0, /C (for instance, since the 
operators are non-local, in all the estimates one should take care of the contribution coming from 
infinity). We plan to address this issue in a future work. 

Assume that ■0 : M" — > M"^ is a smooth globally Lipschitz function, b G a vector, r > is 
a constant and /C a (smooth) non-local translation-invariant elliptic operator of lower order with 
respect to (—A)*, i.e., there exists k E (0, 1) such that 

[/Cy.]c«jM") < ||(-A)VIIl-(M") e C^{R^). 

We consider u : [0, T] x M" — > M a (continuous) viscosity solution to the obstacle problem 

J mm{-ut + ru + b-Vu- {-AYu- ICu,u-'>p} = on(0,r]x]R", 

I n(0) = 0, ^^-'^ 

When s > 1/2, existence and uniqueness of such a solution follows again by standard results on 
obstacle problems. 

Let us now analyze the properties of solutions, as we did before for (|2.2p . 

• Basic properties. We proceed as in Section [3j First of all, as in Lemma 13.11 one can 
approximate solutions to ()5.ip using a penalization method. In this way, all the results of 
Section [3] still hold true: 

- n(t, •) is globally Lipschitz (see Lemma [3.2p : 

- n(t, •) is uniformly semiconvex (see Lemma l3.2p : 

- Ut is globally bounded (see Lemma [3.3p : 

- {—AYu + fCu — ru — b ■ Vu is globally bounded (see Lemma [3.3p . 

In particular, by elliptic regularity for the fractional Laplacian, the L°°-bound on {—AYu + 
K,u{t) — ru{t) — b ■ Vu{t) gives 

- nGL-([0,r],Cf;-°''(M")). 

Hence, since s > 1/2 and ICu is of order < 2s — n, there exists 7 = 7(5, k) > such that 

R ■= -ICu + ru + b-Vue L°°([0, T],Cl^{W')). (5.2) 

• -decay for {—A)^u{t). In this setting, we have: 
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- {-Ayu{t) - R{t) < inside the open set {u{t) > tf;}; 

- {-AYu{t) - R{t) > a.e. on {u{t) = i;}. 

Now, fixed t > and given a free boundary point xq G d{u{t) = ip}, we consider the La-harmonic 
function v{x,y) := u{t,x,y) — -^Y^r^y^~'^, where u{t,x,y) is the La-harmonic extension of u{t). 
Moreover, as in Subsection 14.11 we consider the function 

vix,y) := v{x,y) - tpix). 

Since ut is globally bounded and v{-,0) is semiconvex, all estimates (B1)-(B2) and (B4)-(B5) of 
Subsection 14.11 still hold true, while (B3) becomes 

(B3') limy_j.o+ y"'Vy{x,y) < \R{t,x) — R{t,XQ)\ < C\x — xqP for a.e. x G A, 
limy_^Q+ y"-Vy{x, y) > —C\x — xo\'^ for x G \ A. 

While the proof of Lemma [4. 41 works with no modifications under these assumptions, for the proof 
of Proposition 14.31 we remark that now we do not have limy^Q+ y"'Vy{x,y) > for x G M" \ A, 
which was used to apply Hopf 's Lemma. To overcome this difficulty, in the induction step from 
ko to /cq + 1 one should replace v{x, y) with v{x^ y) + l^'^'^^^^^^^t^it'^Q)) (4~^o^7jyi--a^ and the rest of 
the proof should go through with minor modifications. 
Hence, one still gets sup^^^^.^^ \u{t) — ^| < Cr""*"^* and 

[{-/\yu{t) - R{t)\xu(t)=i, = [{-^yu{t) + /Cn(t) - ru{t) - b ■ Vu{t)]xuit)=^ G Qoc(I^"), 
for some universal exponent a G (0,7). 

• Monotonicity formula and optimal spatial regularity. As in Subsection 14.21 one would like 
to apply a monotonicity formula. However, first of all one has to do a preliminary step: using 
the equation 

Ut + i-AYu = [{-AYu + R]xu=^ -Re L°°{[0, T]; CS,(M")) 

one deduces (thanks to a local variant of (|A21)-(1AH)) that (-A)"n G L'^{[0,T];C^-"^ {M.'')). 
So, by elliptic regularity for the fractional Laplacian, u G Cj"^^'*"''^ (M"), which gives 

i?GL~([0,r],Ci"„+'^/'(M")). 

This allows considering i? as a lower order perturbation when applying the monotonicity formula. 

Now, to apply the monotonicity formula around a free boundary point xq G d{u{t) = ip}, 
one should consider the function w : M" x M+ x R obtained by solving the Dirichlet problem 

L^aU! = 0, 

w{x,0) = [{-Ayu{t,x,0) - R{t,xo)]x{u{t)=i,}(x). 

Since R{t) G C^"'^'^ (W^), we have w > — C'r"^'^/^ on Br{xQ). So, by the monotonicity formula 
one gets [(— A)''n + R{t,x)]xu{t)=tp ^ C'|^c(^")' f^a '•— ™iii{a + 6'^,1 — s}. Then, one can iterate 
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the above strategy, first using the parabolic regularity of dt + {—Ay and the the elliptic regularity 
of (-A)"* to show that R G L°°([0, T]; C(^^+^(M"), and then applying again the monotonicity 
formula. In this way, after finitely many iterations we get 

[{-AYu + Ku-ru-h- Vu]xu=^ G L^{% T], R € L-([0, T], Ci^„7+^(]R")). 

• Parabolic regularity and conclusion. Using Lemma 14.121 the argument in Subsection 14.31 
applied to 

dtu + (-A)^n = [{-Afu + R]xu=^ - R 

allows to extend the regularity result in Theorem 12. II fat least locally in space-time) to solutions 
of ([5TD . 

A Regularity properties of the operator dt + (— A)^ 

In this appendix we describe some important properties of the parabolic operator dt + (—A)*. 

Let us first recall that fractional Laplacian works nicely in Holder spaces: in / S C°(M"') 
then (-A)-*/ € C°+2^(M"), see for instance [El Subsection 2.1]. 

Analogously, the operator dt + (—A)* works nicely in space-time Holder spaces: if vt + 
(— A)*t; = / on [0, r] x M" and z;(0) is smooth (for our purposes, we can assume v{Q) € C^(M"), 
globally Lipschitz, and ||L'^'y(0)||ioo(iRn) + ||( — A)'''y(0)||(-,i-sj'jg„^ < +oo), by classical results on 
multipliers on Holder spaces (see for instance \\?>\ Theorem 2.3] and the proof of |13| Theorem 
3.1]) we get 

"^*"c7,"i''{(0,T]xR") + {(0,T]xR") ~ ^ + Wf^C^'i {{0,T]xW^) ^ (0' 1) (^'l) 

However, for our purposes, we also need to have some regularity estimates when / is only bounded 
in time (but Holder in space). 

Let us observe that we can write the solution in terms of the fundamental solution Tgit^^x) 
of the fractional heat equation. More precisely, using Duhamel formula, we have 

-{-Ayv{t,x) = vt{t,x) - f{t,x) 

= -r,(t) * (-A)MO) + r / dtT,{t-T,x-y)[f{T,y)- f{T,x)]dydT. ^^'^^ 

Jo JR" 

We now claim that the following estimates hold (the proof of them is postponed to Subsection 
[AT] below): 

[ [ dtTs{t - T,x - y)[f{T,y) - f{T,x)]dydT 

Jo JR" 

(A.3) 



< 



I3/{2b),I3-0+ 



■^IIl°°{(0,T];C^{R")) 
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if /5 < 2s, 

dtTsit-T,x-y)[f{T,y)- f{T,x)]dydT - n-/ iiL-ao Tl-c^m")) 

(A.4) 

if /3 > 2s. 

Combining the above estimates with (jA.2p . and using that — rs(t) * (— A)*u(0) is smooth in 
space-time for t > 0, we get 

||(-A)^z;||^,/,,.,,_o+^^^_^j^^„^ < 1 + ll/ll,.c.((o,T];cf(R")) e (0'2s) (A-5) 

ll(-^)^-llio,Up.c.^-o-((o,T]xM") ^ 1 + II/IIlo.((o,t];C.^(m^)) V/5 G [2s, 1) (A.6) 

(At t = the time regularity may degenerate, due to the presence of the term Ts{t) * (— A)*u(0).) 
In particular, using that Vt = f — (— A)*u, we obtain 

"^*"L-((0,T];Cr"+(M")) + "(~^)'^"L-((0,T];Cf-°+(R")) ~ ^ + ll-^llL-((0,T];C;f (R")) ^ ^) 

(A.7) 

A.l Proof of (ig and (1X1) 

Let us recall that t € [0, T], with T < +00. 

To prove ()A.3P and (|A.4p . we use that the fundamental solution rs(l, behaves like ^qq^jsr^j) 
which by scaling implies 

\dttrsit,y)\ < 77^1+2^^——' |V,5tr,(t,y)| < „^,. V , (A-9) 

We will also make use of the following two basics estimates: 
- There exists a constant C > such that, for all h G (0, 1], 



JR" /i^^ + z '^+2'^ ^ 



1/5 

I ) 

There exists a constant C > such that, for all /i > 0, 



(A.IO) 



/ ^ri dr < Cminj-^,-^ [> . (A.ll) 



The proof of both is pretty simple. For instance, to show (jA.lOp . one splits the integral into 
three parts: 
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/ .J"^' , dz< ! \zr-'^ dz<i+ h^/('^)-\ 

To prove (|A.lip . we observe that the bound is trivial if /i > 1, since jj+^s < un+is 

(recall that |t — r| < T < 1). On the other hand, if /i € (0, 1], it suffices to split the integral over 
[0,t — h?^] and \t — h?^ ,t\, and argue as above. 

• Proof of ()A.3p . Let us observe that, for u < t, 

\{-Ayv{t,x)-{-Ayv{T,x)\ 

t f 

dtTs{t -T,x- y)[f{T,y) - f{T,x)]dydT 

dtVsiu -T,x- y)[f{T,y) - /(r,x)] dy dr 

<l I \dtTs{t-T,x-y)\\f{T,y)-f{T,x)\dydT 

+ / dtVsiu - r,x - y) - dtVsit - T,x - y) mm{\x - yf,l}dy dr 
Jo JR" 

= (Tl) + {T2). 

Now, using (jA.Sp . by (jA.lOp applied with h = {t — t) we get that (Tl) is bounded by 

min{|x — y\^ , 1} 



'0 




\n+2s 



dydT< / (l + (t-T)'^/(2s)-n^r 



lu JR" (t — r) 2s -\-\x — y\ 

<{t-u) + it-u 

Concerning (T2), thanks to (|A.8p . (|A.9p . and (jA.lOp with h = (t — t), we can control it by 

1 1 , n+2s ^ , — — min{|x - yf, 1} dy dr 



min 

n 


r 


- u 


[u 


— T 


min "I 


t - 


u 


u — 


r 


-{t-u) 


t - 






u — 


-A 



/ , . I - T I , min{|x- 

[U — t) 2s -\- \x — 

fu-(t-u) ^ _ / 

< 



{l + {u- r)^/(2s)-i^ (l + (n - r)^/(2.)-i^ 



<(t-u) + {t- u)\\og{t -u)\ + {t- uf/^^'\ 
which proves the time regularity of {—/S.yv. 
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Proof of (|A.4p . The proof of the spatial regularity is analogous: we write 
{-Arv{t,x)-{-Ayv{t,z) 

dtVsit - T, X - y)[/(r, y) - f{T, x)] - dtT^it -t,z- y)[f{T, y) - f{T, z)]) dy dr. 

Then, we split the spatial integral over two sets: the region where | |x — z| < \x — 2/|/2} and the 
region where ||x — z| > |x — y|/2|. 

On the first set, since ji+^j^^ and n+is^ are comparable, we can 




+2s 



n+2s 



estimate the integrand by 

\f{r, y) - f{r, x)\\dtrs{t -T,x-y)- dtT^it -T,z-y)\ + |/(t, x) - f{T, z)\\dtrs{t -T,x-y) 
which thanks to (|A.8P and ()A.9P can be bounded by 

\y -X 



y^\X — Z\ 



1 



+ \x — z\ 



(t — r) 2s + |x — y\ 



n+2s 



\x-y\{t-T)'^ + \x-y\^+'^' 
So, using (jA.lip with h = \x — y\ we get 

dtVs{t -T,x- y)[f{T,y) - /(r,x)] - dtTs{t -t,z- y)[f{T,y) - /(r,z)] dydr 




J {\x~z\<\x~y\/2} 



< 



< 



{\x-z\<\x-y\/2} \ \^- ^ 



\X — Z\ 



\x — z\^ min 



+ \x — z\ 



1 



(t - r)^^ + |x - 



dr dy 



l{\x-z\<\x~y\/2} 

< \x — z[ 



1 



{\x-z\<\x-y\/2<l} 



X — yl" ' \x — y\"-+'^^ 
x - y\~"- + \x- zf j 



dy 



{\x~y\/2>l} 



\x - y\ 



~n~2s 



< \x — zl^llog |x — z| I < \x — z\ 



Ve > 0. 



Concerning the integral over the second set, we simply use (jA.SP to bound the integrand by 

\x — y\^ ^ \x — z\^ 

{t - + \x- y|"+2s (t - + \x- z|"+25 

and observing that {\x — z\ > \x — yj/2} C B'^^,j._^^{x) PI B^^^_^^{z) we get 

'dtTs{t-T,x-y)[f{T,y)-f{T,x)]-dtTs{t-T,z-y)[f{T,y)-f{T,z)] 
\x — y\^ 




J {\x~z\>\x~y\/2} 



dy dr 



< 




^B3|:.~.|W {t - t)^ + \x- y|'^+2« 
* 1 



dy dr 



< 



'■B3|x-z|(a;) 



\x - y\ 



n+2. 



(t - r) 2s + \x — y\ 



n+2s 



dr dy 



yf-''dy<\x-zf, 
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where for the last but one inequahty we used again (|A.lip with h = \x — y\. This concludes the 
proof of ()A.4p . 
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